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Jet rates in electron-positron annihilation at 0(af) in QCD 
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We compute production rates for two, three, four and five jets in electron-positron annihilation 
at the third order in the QCD coupling constant. At this order, three-jet production is described to 
next-to-next-to-leading order (NNLO) in perturbation theory while the two-jet rate is obtained at 
next-to-next-to-next-to- leading order (N 3 LO). Our results yield an improved perturbative descrip- 
tion of the dependence of jet multiplicity on the jet resolution parameter, y cu t, particularly at small 
values of y C ut- 
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Jet observables in electron-positron annihilation play 
an outstanding role in studying the dynamics of the 
strong interactions, described by the theory of quantum 
chromodynamics (QCD, l]). The initial experimental 
observation of three-jet events at PETRA in agree- 
ment with the theoretical prediction Q, provided first 
evidence for the gluon, and thus strong initial support for 
the correctness of QCD. Subsequently the three-jet rate 
and related event shape observables were used for the 
precise determination of the QCD coupling constant a s 
(see d, Q for a review), and four-jet observables helped 
substantially to confirm the gauge group structure of 
QCD by firmly establishing the gluon self-coupling [||. 

Jets are defined using a jet algorithm, which describes 
how to recombine the momenta of all hadrons in an event 
to form the jets. A jet algorithm consists of two ingre- 
dients: a distance measure and a recombination proce- 
dure. The distance measure is computed for each pair 
of momenta to select the pair with the smallest separa- 
tion. This pair of momenta then is combined according 
to the recombination procedure into a joint momentum, 
if its separation is below a pre-defined resolutionparame- 
ter j/cut- Improving upon the JADE algorithm 7], which 
uses the pair invariant mass as distance measure, several 
jet algorithms have been proposed for e + e _ collisions: 
Durham Geneva P and Cambridge [icj |. Among 
those, the Durham algorithm has been the most wi dely 
used by experiments at LEP [H OS 0, and SLD pf, 
as well as in the reanalysis of earlier data at lower ener- 
gies from JADE [jj]. 

The Durham jet algorithm clusters particles into jets 
by computing the distance measure 



Vij,D 
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for each pair of particles, E v - ls denotes the energy 

sum of all particles in the final state. The pair with the 
lowest yij,D is replaced by a pseudo-particle whose four- 
momentum is given by the sum of the four-momenta of 
particles i and j ('E' recombination scheme). This pro- 
cedure is repeated as long as pairs with invariant mass 
below the predefined resolution parameter yij.D < J/cut 



are found. Once the clustering is terminated, the remain- 
ing (pseudo-)particles are the jets. It is evident that a 
large value of y cu t will ultimately result in the clustering 
all particles into only two jets, while higher jet multi- 
plicities will become more and more frequent as y cu t is 
lowered. In experimental jet measurements, one there- 
fore studies the jet rates (jet cross sections normalized 
to the total hadronic cross section) as function of the jet 
resolution parameter y cu t- 

The theoretical prediction of jet cross sections is made 
within perturbative QCD, where the same jet algorithm 
is applied to the momenta of final state partons. The 
QCD description of jet production is either based on a 
fixed-order calculation, which uses exact parton-level ma- 
trix elements (including higher order corrections if avail- 
able) for a given jet multiplicity, or by a parton shower, 
which is based on the leading-order matrix element for 
two-jet production only, and generates higher multiplic- 
ities in an iterative manner, thereby accounting only for 
the leading logarithmic terms from parton-level processes 
with higher multiplicity. Depending on the jet multi- 
plicity, higher perturbative orders correspond to different 
powers of the QCD coupling constant: the leading order 
prediction for n-jet production is proportional to a"~ 2 . 
So far, fixed-order calculations were available up to next- 
to-next-to-leading order (NNLO) for two jets fl^.l^. [l9l . 
up to next-to-leadin g or der (NLO) for three [20 , |2U |22| 
and four jets [23l l24l l25l |2(| . For five and more jets, only 
leading order calculations were available [13, Hs HH ■ F° r 
jets involving massive quarks, NLO results are available 
for three-jet final states [3fj| . 

Calculations based on parton showers, incorporated 
in multi-purpose event generator programs [3l|, [32|, [H| , 
provide a satisfactory description of multi-jet production 
rates. Since these programs contain many tunable phe- 
nomenological parameters, their predictive power is how- 
ever very limited. 

In this letter, we present the first calculation of NNLO 
corrections to three-jet production and the next-to-next- 
to-next-to-leading order (N 3 LO) corrections to two-jet 
production in e + e~ annihilation. Together with the pre- 
viously available NLO corrections to four-jet production 
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FIG. 1: Perturbative fixed-order description of the three-jet 
rate at Q = Mz, compared to data obtained with the ALEPH 
experiment [llj. Experimental errors are too small to be vis- 
ible on the figure. 



and the leading-order description of five-jet final states, 
these are used for a fully consistent perturbative descrip- 
tion of e + e~ — > jets at order a^. in perturbative QCD. 

The calculation of the ai? corrections for three-jet pro- 
duction is carried out using the newly developed parton- 
level event generator program EERAD3 which contains the 
relevant matrix elements with up to five external par- 
tons [13, [H, [H, H(| . Besides explicit infrared divergences 
from the loop integrals, the four-parton and five-parton 
contributions yield infrared divergent contributions if one 
or two of the final state partons become collinear or soft. 
In order to extract these infrared divergences and com- 
bine them with the virtual corrections, the antenna sub- 
traction method [HI, H3] was extended to NNLO level (38j 
and implemented for e + e~ — > 3jets and related event- 
shape variables into EERAD3. The analytical cancel- 
lation of all infrared divergences serves as a very strong 
check on the implementation. 

Initial results obtained with EERAD3 on NNLO correc- 
tions to event shape observables were reported in (4^| and 
applied in the extraction of the strong coupling constant 
from LEP data in [4l[. Since the program provides the 
full kincmatical information for each event, it can also 
be used to simultaneously compute the production cross 
sections for three, four and five jets through to 0{a? s ) 
for any infrared-safe jet algorithm and as function of the 
jet resolution parameter. The jet rates are then defined 
by normalizing the multi-jet cross sections to the total 
hadronic cross section computed at the same order. 

The four-jet [H H, [H, ^ and five-jet rates [13 were 
known previously to O{o i s ). Our major new result is 
the three-jet rate to this order, which corresponds to 
NNLO in the perturbative expansion. Figure [1] displays 
the three-jet rate at LEP1 energy Q = Mz as function 
of the jet resolution y cut at LO, NLO, NNLO. At NNLO, 
the denominator has been expanded, as described in [40j | 



to contain only terms up to O(a^) in the jet rate. The 
theoretical uncertainty band is defined by varying the 
renormalization scale fj, in the coupling constant in the in- 
terval Mz /2 < n < 2 Mz, and the world average value Q 
a s (Mz) = 0.1189 is used, consistently evolved to other 
scales at each order. The fixed-order theoretical predic- 
tions for three-jet rate become negative for small values of 
?/ cu t, where fixed order perturbation theory is not applica- 
ble due to the emergence of large logarithmic corrections 
at all orders, requiring resummation [1,[43|. We therefore 
restrict our comparison to y cu t > 10 -4 , although data at 
lower jet resolution parameters are available. 

For large values of y C ut, Vcut > 10~ 2 , the NNLO cor- 
rections turn out to be very small, while they become 
substantial for medium and low values of y cu t- The maxi- 
mum of the jet rate is shifted towards higher values of y cu t 
compared to NLO, and is in better agreement with the 
experimental observation. The theoretical uncertainty is 
lowered considerably compared to NLO. Especially in the 
region 1CP 1 > y cu t > 1CP 2 , which is relevant for preci- 
sion phenomenology, one observes a reduction by almost 
a factor three, down to below two per cent relative uncer- 
tainty. Since the error band in this region is barely visible 
in the plot, we display the relative theoretical uncertainty 

_ max A1 (cr(/i)) - min M (cr(/i)) 
2<t(h = M z ) 

at NLO and NNLO as an inset. The relative uncertainty 
on the LO calculation is constant at 10.2%. 

The fixed-order NNLO description is still above the 
data at low jet resolution, where the convergence of the 
perturbative series is spoilt by large logarithms of y cut at 
all orders, and where a resummation should be carried 
out Q. Furthermore, the theoretical parton-level predic- 
tion is compared to hadron-level data, thereby neglect- 
ing hadronization corrections, which may also account 
for part of the discrepancy. 

To compute the jet rates with different multiplicities, 
it is more appropriate to normalize all jet cross sections 
to the total hadronic cross section corrected to third or- 
der [H|] in the QCD coupling constant, Q(a^). We con- 
sistently neglect numerically small QCD singlet contri- 
butions at this order, which were found to contribute at 
most one per cent [13j to the total coefficient of the 0( a s) 
correction, and which are equally small in the individual 
jet multiplicities [H|. The total hadronic cross section 
is made up from the sum over all jet multiplicities. At 
0{oP 8 ), this sum runs from two-jet through to five-jet fi- 
nal states, such that the corresponding jet rates must add 
to unity. Consequently, our calculation yields the N 3 LO 
expression for e + e~ — > 2 jets as a by-product. It is in- 
teresting to note that some earlier NNLO calculations of 
the two-jet rate [l8l [l^ | were essentially exploiting the 
same feature at 0{a 2 s ). 

Figure [2] shows the parton-level theoretical predictions 
for the jet fractions at first, second and third order in 
the strong coupling constant, compared to experimental 
hadron-level data from ALEPH [Tlj] . 
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FIG. 2: Jet rates at first, second and third order in the strong coupling constant, compared to data from ALEPH [III ]. The 
rates are normalized to the total hadronic cross section at that order. 



By comparing the three plots, we observe that there 
is systematically improved agreement for each of the jet 
rates as the order of perturbation theory increases. At 
each order a new multi-jet channel opens up, e.g. the 
five-jet rate at 0(al), which is positive definite and es- 
sentially monotonically increasing at small y cut . Since all 
jet rates are normalized to unity, the new five-jet channel 
has the effect of reducing the contribution to the two-jet, 
three-jet and four-jet rates, in the region of log 10 (j/ cllt ) 
where the five-jet rate contributes. One very clear ef- 
fect is to cause the turnover in the four-jet rate (which 
is not present at 0{a 2 s )). A second effect is to add more 
structure to the shape of the two- and three-jet rates, 
which lie much closer to the data for log 10 (?/cut) < —2.5. 
Of course, the effect of the higher order corrections also 
extends to larger values of y cu t , due to the different con- 
tributions of the two-loop virtual and virtual-radiation 
graphs to the three- and four-jet rates, as well as the way 
that the double radiation contribution interacts with the 
jet algorithm, and through the normalization to the total 
hadronic cross section. This is visibly less dramatic, but 
by adding more structure to the theoretical prediction, 
enables a better description of the data. 

Previous experimental studies of multi-jet production 
rates compared only with standard leading-order parton 
shower event generator programs, which yielded a good 
description of the data at the expense of large hadroniza- 
tion corrections 0, [ll| . In the light of our new results, 
this issue should be carefully reexamined within fixed- 
order perturbation theory. 

In this letter, we reported on the NNLO QCD correc- 



tions to the three-jet production rate at parton-level in 
e + e~ annihilation, which is the first genuine NNLO cal- 
culation of a jet production rate at particle colliders. We 
observed that (hadron-level) experimental three-jet data 
are described considerably better in shape and normal- 
ization, and over a wider range in y cu t, than at NLO. 

At the same order in the strong coupling constant, a^, 
we describe four-jet production at NLO and five-jet pro- 
duction at LO, reproducing earlier results. By combining 
those and normalizing to the total hadronic cross section 
at this order, we obtained the two-jet rate to N 3 LO in 
perturbation theory as a by-product. We observe that 
with increasing order in the strong coupling constant, 
the multi-jet rates are better described over an increas- 
ing range of resolution parameters. Our results clearly 
highlight how perturbative QCD successfully describes 
jet production rates at the parton-level. 
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